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ABSTRACT 

Wc compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) fre- 
quency spectrum of a polar mountain produced b y magnetic buria l on an accreting 
neutron star. Applying the formahsm developed bv lHellsten fc Spiei (jlOTOD . extended 
to include gravity, we solve the singular eigenvalue problem subject to line-tying 
boundary conditions. This spectrum divides into an Alfven part and a cusp part. 
The eigenfunctions are chirped and anharmonic with an exponential envelope, and 
the eigenfrequencies cover the whole spectrum above a minimum uj\ow For equilibria 
with accreted mass 1.2 x 10~^ < M^/Mq < 1.7 x 10~^ and surface magnetic fields 
10"'^"'^ < i3*/G < 10^'^, wiow is approximately independent of B*, and increases with Ma- 
The results are consistent with the Alfven spectrum excited in numerical simulations 
with the ZEUS-MP solver. The spectrum is modified substantially by the Coriolis force 
in neutron stars spinning faster than ~ 100 Hz. The implications for gravitational 
wave searches for low-mass X-ray binaries are considered briefly. 

Key words: accretion, accretion disks - stars: magnetic fields ~ stars: neutron - 
pulsars: general 



1 INTRODUCTION 

Radio and X-ray observations suggest that the magnetic 
dipole moments, jj., of neutron sta rs in accreting binaries 
decre a se with accreted mass, Mg (|Taam fc van de Heuvell 



1 19861 : Ivan den Heuvel fc Bitzarakil 19951 '). One physical 



mechanism capable of reduci ng ^ by the amount observed is 
magnetic screening o r burial ('Bisnovatvi-Kogan fc Rom berg " 
1974 iRomanil I1990I: IXonar fc Bhattacharva 1997 : Zhang 



1998: Melatos fc Phinnev 2001: Choudhuri fc KonaJ 1200: 



Payne fc MelatosI |2004| : iLovelace et al.l 120051 ). Magnetic 
burial occurs when accreting plasma, flowing inside the 
Alfven radius, is channelled onto the magnetic poles of 
the neutron star. The hydrostatic pressure at the base 
of the accreted column overcomes the magnetic ten- 
sion and the column spre ads equatorwards, distort ing the 
frozen- in mag netic flux (|Melatos fc PhinnevI 1200 ll ). Self- 
consistent magnetohydrodynamic (MHD) equilibria respect- 
ing th^_^iix;^eezing, ideal- MHD constraint were computed 
by iPavne fc Melatoa (|2004l ) , who found that the magnetic 
field is compressed into an equatorial belt, which confines 
the accreted mountain at the poles. A key result is that 
/i drops significantly once Ma exceeds the critical mass 
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Mc ~ 10~ Mq. This characteristic value exceeds simple es- 
timates based on local MHD force balanc e at the polar cap , 
i.e. without the equatori al magnetic belt l|Brown fc Bildster] 
1 19981 : iLitwin et al1l200ll ) . 

Counterintuitively, magnetic mountain equilibria 
prove to be marginally stable (jPavne fc Melatoa l2007l : 
I Vigelius fc MelatosI 120083 ). An axisymmetric mountain is 
susceptible to the undular submode of the Parker instabil- 
ity, but the instability is transitory, saturating after ~ 10 
Alfven times to give a nearly axisymmetric, mountain-like 
state, which oscillates in a superposi tion of small-amplitude , 
global, Alfven and acoust ic modes (jPavne fc Melato3l2007l : 
IVigeUus fc MelatosI [20083 ). The magnetic line-tying bound- 
ary condition at the stellar surface is crucial in stabilizing 
the mountain. 

The Alfven and acoustic oscillations help to shape 
the gravitati onal wave spectrum emitt ed by low-mass X- 
ray binaries (|Vigelius fc MelatosI l2008al ) . furnishing a new 
observational probe of the surface magnetic structure of 
neutron stars. The oscillations may also manifest them- 
selves as small, Hz-to-kHz variations in the X-ray pulse 
shape. Accordingly, a linear eigenmode analysis of the ideal- 
MHD spectrum is required to take full advantage of fu- 
ture gravitational- wave and X-ray timing experiments. Here, 
we take a first step by computing the continuous part 
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of the ideal-MHD spectrum analytically for axisymmet- 
ric magnetic mountain s, closely following the approach of 
iHellsten fc Spied (|l979[ l. 

Realistically, however, the detection of the 
gravitational- wave imprint from mountain oscillations 
will not be possible in the near future. While Advanced 
Laser Interferometer Gravitational Wave Observatory 

i Advanced LIGO) may detect the unperturb ed mountain 
Watts et all l2008l : IVigehus fc Melatoj [ioOSal ). the detec- 
tion of modulations of the main signal from mountain 
oscillations will have to await next-generation interferome- 
ters. On the other hand, X-ray-profile changes c an currently 
be measured with an accur acy of ~ 1 per cent (|Muno et al.l 
|2002| : iHartman et al.ir2008l ) and we know that the accretion 
rate changes by ~ 10 per cent per day. Unfortunately, the 
X-ray fluctuations will be only poorly frequency-matched 
to the mountain oscillations in general. 

Generally, an inhomogenuous MHD c o nfigu- 
ration supports linear eigenmodes (|Lifschit j 1 19891 : 
iGoedbloed fc Poedtd |2004 ). whose frequency spectrum 
divides into a discrete and a continuous part. Discrete 
eigenvalues are fixed by the boundary conditions. Contin- 
uous eigenvalues arise from singularity in the underlying 
Sturm-Liouville problem, which allows the boundary 
conditions to be satisfied for any eigenvalue within a range. 

In earlier work, stochastically excited mountain oscil- 
lations were investigated numerically by perturbing equi- 
librium s tars with different Ma in the ideal-MHD solver 
ZEUS-MP jHaves et al.ll2006h. and computing the sp ectrum 
jPavne fc Melatod l2007l : IVigehus fc MelatosI boOScT ) . How- 



ever, this approach is restricted to the subset of the full 
MHD spectrum resolved by zeus-mp and is computation- 
ally expensive. In this article, we attack the problem ana- 
lytically. The article is organised as follows. In section [21 we 
introduce curvilinear field line coordinates to describe the 
equilibrium and establish the associated metric. We derive 
the linearized, ideal-MHD equation s in these coordinates in 
section [3] extending the analysis bv IHellsten fc Spies! (|l979l ) 
to include the gravitational field of a central point mass. 
The continuous frequency spectrum and the corresponding 
eigenfunctions are evaluated in section |4j as a function of 
Ma and the magnetic field strength before burial. We con- 
clude by discussing the implications for gravitational wave 
observations of accreting millisecond pulsars in section [5] 



2 HYDROMAGNETIC EQUILIBRIUM 

The equilibrium structure of a magnetically confined moun- 
tain in ideal MHD is described by the force balance equation 



VP - (V X B) X B = Ng, 



(1) 



supplemented by V ■ B = and an equation of state P{N), 
which we take to be isothermal; P = c^A''. Here, P, B, N, 
g, and Cs denote the pressure, magnetic field, mass density, 
gravitational acceleration, and isothermal sound speed, re- 
spectively. 

If we introduce a cylindrical coordinate system (r, ip, z) 
and assume axisymmetry, we can write 



where ^ is the magn etic flux, measured in G cm . 
iPavne fc Melatod (|2004l ) computed unique, self-consistent, 
MHD equilibria by solving ^ in spherical coordinates sub- 
ject to the flux- freezing constraint of ideal MHD. In this 
article, we convert these solutions to cylindrical coordinates 
before constructing the associated field line coordinates. 

Throughout this article, we work in cgs-like units, such 
that ^0 = 1- Furthermore, we normalize the isothermal 
sound speed and the gravitational constant, viz. Ca = G = 1. 



2.1 Field line coordinates 

A magnetic mountain is created by continuously deforming a 
dipole magnetic field during the accretion process. The flux 
surfaces are closed at all times, a topology that is preserved 
even for a uniform ^ 0. We can therefore introduce or- 
thogonal curvilinear coordinates {^,rj,ip), called field line 
coordinates, which "follow" the shape of the flux surface. In 
these coordinates, 77 measures arc length along a magnetic 
field line, normalized to the domain ^ 77 ^ 1, where 7; = 
corresponds to the footpoint at the surface, and rj — 1 corre- 
sponds to the outer radial boundary or, when the field line 
is closed, to the surface, 'i' is the fiux function in Q and ip is 
the usual azimuthal angle in cylindrical co ordinates. In the 
field line coordinates, the metric becomes (|Hellsten fc Spie3 

ds'' = — i-d^'^-f J'bW+'-'V, (3) 
with the Jacobian J given by 

J={B-Vri)-'. (4) 

Note th at, for B^ = 0, we have J = B^^ . However, in keep- 
ing with lHellsten fc Spied (| 19791 ) we uphold the more general 
notation in order to facilitate the inclusion of a toroidal field 
component in future work. 

The contravariant components of the force balance 
equation ([T} in the 9 and rj directions read respectively 



dP i d , ^^2-. ,r 

and 

— = AiN. 
or] 



(5) 



(6) 



The gravitational acceleration is directed radially inward to 
the centre of the star (mass M). Ignoring self gravity, we 
have 



g = AoBis, + Aie,,, 
with 

. ... . dr dz\ 

Ao = g ■ = --r^, — t^^tttt; ( r— + z- — 1 



M / 



and 



M 



dr dz 
r— — h z — 



(7) 
(8) 

(9) 



\/(p X V*, 



(2) 



(r^ + z^)^/^ \ dr] dr] 

IPavne fc Melatod (|2004l ) took g to be constant to sim- 
plify the analysis, but we prefer to use Eqs. (|8]) and ([U 
to allow direct c o mpar ison with the numerical results of 
iPavne fc Melatod (|2007l ) and I Vigelius fc Melatod (|2008d ). 
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For later comparison, we note that ( O and Q become 
forma lly identical to Eqs. (13) and (14) of lHellsten fc SpiesI 
l|l979l ) for a stationary fluid in the absence of gravity, rotat- 
ing around the z axis with constant angular velocity SI, if 
we substitute 



and 



(10) 



(11) 



Of course, we cannot write the centrifugal force in the form 
0, as it is not radial, so the substitution pup and (|ll|l is 
algebraic, not physical. 



2.2 Accreted magnetic mountain 

Throughout this paper, we study magnetic mountain equi- 
libria on the surface of a curvature-downscaled star with ra- 
dius R'^ = 2.7 X 10^ cm and mass M', = 1.0 x lO'^^M©. The 
downscaling transformation preserves the equilibrium shape 
of the mountain [exactly in the s mall Mg limit and a pprox - 
imately in the large Mg lim it; see lPavne fc MelatosI (|2004 ): 
IVigelius fc Melatod (|2008cl )] , as long as the hydrostatic scale 
height ho = 53.82 cm keeps its original value for a realistic 
star. Base units are Mo = 8.1 x 10^* g, A^o = 5.2 x 10^^ g 
cm"^, Bo = 7.2 X 10^'' G, and ro = 5.4 x 10"'^ s, for mass, 
density, magnetic field, and time respectively. In order to 
upscale the frequencies back to a r ealistic neutr on star, we 
employ the relation cx {ho/R*f (|Pavne fc Mclatos,200^ . 
While, strictly speaking, this relation only applies to waves 
travelling latitudinally, we note that the field lines are pre- 
dominantly parallel to the neutron star surface (Fig. [T] left 
panel) and the error will be small. 

An equilibrium configuration with Mg = Mc — 1.2 x 
10"'' Mq is displayed in Fig. [T] Mc denotes the critical ac- 
creted mass beyond whi ch the magnetic dipol e moment 
is substantially reduced (jPavne fc Melatod [iooi ) . The left 
panel shows the density contours (dashed) and the magnetic 
field lines (solid), i.e. the flux surfaces *!/ = const, projected 
into a meridonial plane. The equilibrium configuration ex- 
tends over ^ X — {R — i?,)//io 10 and ^ 9 ^ n, where 
R measures the radius in spherical polar coordinates, and 
9 is the colatitude. North-south symmetry is assumed. The 
mountain is confined to the magnetic pole and the distorted 
magnetic belt is clearly visible at 9 ^ 0.5 and r ^ 52. The 
right panel shows the same plot in the ^-77 plane. Of course, 
the field lines are projected onto straight lines in this plot. 
The density contours appear distorted since rj is normalized 
to the domain ^ rj ^ 1. 

Contours of B are plotted as dotted curves in both pan- 
els of Fig. [T]for log„)(B/Ba) = -7, -6, -5.5. The magnetic 
belt with its enhanced magnetic field is clearly visible be- 
tween 0.6 ^ ^ ^ 1.4 in the left panel. The dotted curves 
trace out isosurfaces of magnetic pressure and are therefore 
useful for visualizing the Lorentz force confining the moun- 
tain. 



3 GLOBAL LINEAR MHD OSCILLATIONS 

We now consider the behaviour of small-amplitude pertur- 
bations of the magnetic mountain equilibria described in sec- 
tion (2] The linearized equations of ideal MHD are projected 
onto the field line coordinate system in section [STT] The sin- 
gularities in these equations, which determine the form of 
the continuous MHD spectrum, are located in section \X2\ 



3.1 Equations of motion 

Follow ing the notation and approach of iHellsten fc Spiej 
(| 19791 ). we expand the velocity and magnetic field pertur- 
bations in terms of their contravariant vector components: 



V = V + 7J Br. 



b = f)*eip 



(12) 
(13) 



As in any curvilinear coordinate system, the components do 
not have the same units in general, e.g. «* and v'^ have units 
G cm^ s~^ and s~^ respectively, since and have units 
(cm G)~^ and cm respectively. The total pressure perturba- 
tion is defined as 



Tv = p + B- b^n + JB b 



2,v 



(14) 



where p denotes the hydrostatic pressure perturbation. 

We then Fourier decompose the perturbed variables 
with respect to time and </?, e.g. v oc exp[i{—ujt + mip)], re- 
calling that the equilibrium is Eissumed to be axisymmetric. 
Thus we can write down the components of the linearized 
momentum balance equation. 







Lo'^Nv^ \LoAQnx 



9* J a* 



J d-q \r'^B 



= uj^NJ'b'^v'^ - icuAm 
d-K . d 
drj drj 



= ui'^Nr^v'^ — uimn — (r'^b'^) , (17) 

J or] ^ ' 

the components of the linearized induction equation, 
db/dt = V X {v X B), 



i-^"- + =0, 

lujb*^ + ^-H— = 0, 
J or] 



(18) 
(19) 
(20) 



and the linearized mass continuity equation, dn/dt + V 
(Nv) = 0, 

.aw" 



iLon = iLoJNb^ + N' 



drj 



Equations (|15|l - (|2ip are derived by writing down the vector 
operators in the curvilinear coordinate system, e.g. 
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Figure 1. Hydromagnetic structure of a magnetically confined mountain with Ma = Mc in spherical polar and field line coordinates. The 
left panel displays the density contours (dashed) \og-^Q{N/No) = -13, -12, -11, -10.7, -10.5, -10.3, with A^o = 5.2 X 10^^ g cm-^, the 
magnetic field strength contours (dotted) log^Q(-B/-Bo) = —7, —6, —5.5, with Bq = 2.55 X lO'^^ G, and the magnetic field lines (solid) with 
footpoints at R = Rt and 9 = 0.10, 0.12, 0.15, 0.20, 0.39, 0.79, in the R-9 plane, where (R, 6, ip) are standard spherical polar coordinates. 
The right panel presents the same information in the SP-r; plane, where ^ is normalized to the fiux surface that closes at the inner edge 
of the accretion disk ^a, and rj is the arc length along a field line, normalized to ^ ^ 1. 



■.Of 



df 



1 df 



d'f 



dri 



(22) 



3.2 Singular eigenvalue problem 

Equations (I15ll-( [21l) can be cast into the form 
(|Hellsten fc Spied |W79 ) 

C^ + DX + EY=0, (23) 
FY + GX = 0. (24) 

In (|23[) and (|24|) . X is a 2- vector with elements and tt, 
while Y is the five-vector (u'', u"^, fe*, fe'', fo*')''". C and D are 
matrices, while E, F, and G are difi'erential matrix operators 
involving ordinary derivatives in rj. Using (|24p . Y can be 
eliminated from (|23|l to give 



dX 



+ C (D EF G)X 



0. 



(25) 



The continuous part of the frequency spectrum consists 
of the values of uj for which psp becomes singular at a flux 
surface. This happens when either 



CX = 



(26) 



N 







B ■ V 



B'^+P 





B-V]+A 



with 



A = B-V 



WTp 



BAiB ■ ViV + 



B2 + P 



(29) 



(30) 



(31) 



The eigenvalue equation (|28|l decouples into two second- 
order, ordinary differential equations (involving derivatives 
B • V = Bd/dri), one for the Alfven continuum and 
one for the cusp cont inuum. The problem is self-adjoint 
iHellsten fc Spies! Il979l ). so the eigenvalues uj^ are real. In 
the absence of gravity (^4 = 0), the operator L is negative 
and both continua are real {uj'^ ^ 0). Gravity provides a 
constant offset A, such that the Alfven continuum is stable 
while the cusp continuum may be unstable. We note that 
the cusp continuum is the two-dimensional equivalent of the 
slow magnetosonic continuum in a one-dimensional, gravi- 
tating plasma slab (see section IT2l and references therein for 
a full discussion). 



FY = 



(27) 



have nontrivial solutions at some 'i' = 'i'o- The flow contin- 
uum, equation (|26|l . has only the trivial eigenvalue ui = 0. 
Equation (|27|l defines the Alfven and cusp conti nua. It can 
be rewritten in the form (|Hellsten fc Spieal 19791 ) 

Lo'^Pr] + Lr] = 0, (28) 

where r] = {Jv"^ , v'^)'^ and the matrices P and L are given 
by 



4 CONTINUOUS SPECTRUM OF A 
MAGNETIC MOUNTAIN 

4.1 Algorithm 

We solve equation (|28p numerically for the continuous spec- 
trum by following the procedure below. 

(i) We compute the axisymmetric Grad-Shafranov equi- 
librium, employing an iterative relaxation algorithm devel- 
oped bv lPavne fc MelatosI (|2004l ) (see section [2^. 
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(ii) The magnetic field equation, dr/diy = B{r), is 
integrated via a fourth order Runge-Kutta algorithm 
ijPress et al.lll98^ to obtain field lines r = r(r;,\E'o), z — 
z{r], ^o) starting from different footpoints at the stellar 
surface. 

(iii) The field values and their derivatives are evaluated 
along the field lines. The entries in P and L are computed 
by spline interpolation in 77. 

(iv) We find uj^ and the associated eige nfunction by using 
a shooting algorithm (jPress et al.l Il986l ) to integrate l|28|) 
separately for the Alfven and cusp continua. The integrator 
is fourth-order Runge-Kutta with adaptive step size control. 

The boundary conditions for the eigenfunctions along 
any field line (^o,V^o), are: 

(i) Jv'^ir, = 0) = w^(77 = 0) = 0, v^'i'n = 1) = v'^iri = 1) = 

when the field line is closed, or 

(fi) d^{Jv'')iv = 1) = d^v^ir^ = 1) = 

when the field line leaves the integr ation area, consis- 
tent w ith IPavne fc Melatod (|2007l ') and IVigehus fc MelatosI 
l|2008d ). The boundary conditions at the stellar surface en- 
force line tying, while the zero-gradient outer boundary con- 
dition crudely approximates the mag netosphere- accretion 
disk c oupling [cf. the discussion in IVigelius fc MelatosI 
l|2008d )]. 

We present the eigenfunctions Ju'' and v'^ for four se- 
lected fieldlines in section 14.21 two leaving the integration 
volume and two returning to the surface. The continuous 
spectrum and its dependence on the accreted mass Ma and 
surface magnetic field strength B* are discussed in section 
K3\ 



4.2 Eigenfunctions 

We preface this subsection by briefiy reviewing the physical 
origin of the continuous spectrum in a plane-parallel, gravi- 
tating (but not self-gravitating) plasma slab. This analogous 
system can be treated analytically and is helpful when inter- 
preting the result s for a magnetic mountain. H ere, we follow 
the exposition in lGoedbloed fc PoedtsI (|2004l '). 

Consider an infinite slab in Cartesian coordinates 
{x,y,z), whose magnetic flux surfaces are perpendicular to 
the gravitational acceleration (directed along the a;-axis). 
We render the problem one-dimensional by assuming that 
all equilibrium quantities N{x) and B{x) — B{x)ey depend 
only on the height x or, equivalently, on \E'(x). In this case, 
equation (|25|) exhibits a genuine singularity, i.e. the asso- 
ciated eigenfunctions become singular, when the eigenvalue 
uj^ equals the local Alfven or slow magnetosonic frequency, 
i.e. either uj^ = uj^i^) ~ C^o ■ B)^/N or lo'^ — tosi^) ~ 
Nuj\{x)/{N + B^). The horizontal wave vector feo is the 
projection of k onto the magnetic flux surface. All eigen- 
functions can be written as ^{x) — ^o(^) exp[i(fcyj/ -I- kzz)], 
with fco = {ky + k^Y^^. It can be shown that w* is then 
square integrable, involving a logarithmic singularity, while 
is not square integrable for the Alfven continuum, and 
v'" is not square integrable for the slow modes. In addition, 
(|25p exhibits an apparent singularity, when ui^ equals the 
local magnetosonic turning point frequencies, i.e. the slow 



and fast magnetosonic frequencies for k^ = 0. In this case, 
the eigenfunctions remain flnitfEI. 

The mountain equilibrium exhibits logarithmic singu- 
larities in X an d Y f or the Alfven and cusp continuum 
iHellsten fc Spiej (|l979l ). Naturally, in a realistic scenario 
perturbations are damped by nonideal effects such as vis- 
cosity and resistivity and singular eigenfunctions cannot 
arise. For example, resistivity will limit the maximum frac- 
tional amplitude of a magnetic perturbation to SB/B ^ 
ln[(rcr;)^''^], where rc is a characteristic time scale and r] 
is the resistivity. If we choose tq = ta = 2.2 x 10~^ cm 
s~i and 77 = 1.3 X 10^^ (|Vigelius fc MelatosI [2008bh we find 
5B/B < 30. Of course, the linear approximation breaks 
down at such a high fractional amplitude. 

We are now ready to apply these ideas to the contin- 
uous spectrum of a two-dimensional, magnetically confined 
mountain. In a two-dimensional system, all eigenfunctions 
are functions of and 77. Fig. [2] (right panels) shows the 
eigenfunctions v^{ri) for the Alfven continuum [lower row 
of Eq. (|28|) ] of an accreted mountain with Ma = Mc. The 
foot points of the field lines are at ^ = 0.1,0.2,0.6,1 rad, 
corresponding to *o/*a = 0.03,0.12,0.96,2.12. Here, is 
the flux surface that closes at the inner edge of the accretion 
disk. The fleld lines are traced out and labelled in the left 
panel. The right panels show v'^{ri) for the same field lines 
(top to bottom) for different numbers n = 0, 1, 30, where n 
equals the number of nodes in the eigenfunction. 

Field lines CD and ® leave the integration volume, while 
field lines CD and ® close back onto the surface. Field line CD 
exhibits only small curvature. Consequently, u*^ is sinusoidal. 
Its envelope increases exponentially oc exp(1.577) (best seen 
in the rightmost panel). Its wavelength also increases with 
ri, which can be understood by noting that the Alfven speed 
DA oc exp(4.5r7) increases towards the outer boundary. In 
contrast, field line ® runs parallel to the stellar surface until 
it takes a sharp turn at 77 ~ 0.8. Thus, t;*^ is a sine wave 
whose amplitude and wavelength remain constant until the 
curvature term B ■ V changes at the bend, sharply increasing 
the amplitude and the wavelength. 

The closed field lines (D and ® behave like field line CD. 
Line (D shows a spiky feature at 77 = 0.5 (equator) due to 
the change of curvature locally. Line ® is almost a perfect 
sine wave for all 77. Recall that, in a plane-parallel slab, N 
and B (and consequently va) remain constant for all 77. The 
eigenfunctions are hence pure sine waves with the disper- 
sion relation da ~ uj/k. We discuss this dispersion relation 
quantitatively in the next section. 

The eigenfunctions Jv^ of the cusp continuum, dis- 
played in Fig. [3] (right panels), behave similarly. Field lines 
CD and (D run parallel to the stellar surface until they bend 
radially outwards. Along these field lines, Jd'' is essentially 

^ An even simpler system is obtained by considering an expo- 
nentially stratified atmosphere with uniform sound and Alfven 
speeds, i.e. constant o;^ and tjg. Under these circumstances, the 
continuous spectra degenerate into a single point each, which 
are cluster points of the discrete spectra. This case includes the 
Parker instability l|ParkeJll967l : lMouschovia3ll974) . The Parker 
instability is responsible for the transient, three-dimensional, 
ideal-MHD relaxation of an initially axisymmetric, magnetically 
confined mountain observ ed in previous numerical simulations 
llVieelius fc MelatodlioOSd) . 



6 M. Vigelius et al. 




n=0 n=1 n=30 



Figure 2. Toroidal velocity eigenfunctions of the Alfven continuum for a magnetic mountain with accreted mass Ma = Mc = 
1.2 X 10~*Mq. The foot points of the field lines (D, ®, ®, and ® are at 6» = 0.1, 0.2, 0.6, 1 rad (left panel, tracing the field lines of the 
mountain in the meridional plane). The right panels graph ii*^ against arc length, r], along the four field lines (top to bottom) for node 
numbers n = (left), 1 (middle), 30 (right), ho = 53.82 cm is the characteristic scale height (units of axes in left panel) and Cs = 10* 
cm s~^ is the sound speed, rj is normalized to the length of the field line, such that ^ ^ 1. 



sinusoidal, with a sharp rise in amplitude at ~ 0.4 and 
77 ~ 0.8 respectively. This behaviour can be explained by the 
presence of the coefflcient NB^ /(B^ + N), which contains 
the partial hydrodynamic pressure (the factor B^ conies in 
since we compute Jv^ instead of v^), which is absent in the 
equations for the Alfven continuum. A comparison with Fig. 
[1] shows that B and TV flatten out considerably at 77 = 0.4 
{rj = 0.8) for field line ® ((D). Consequently, the derivative 
of the coefficient is steep up to that point and becomes neg- 
ligible thereafter. Field lines (D and ® behave essentially like 
their Alfven counterparts, showing the same sharp spike at 
»7 — 0.5 as lines CD and ®. Remember that the field lines ® 
and ® cover only one hemisphere, i.e. r] — 1 corresponds to 
the equator (see Fig. [T|. 



We summarize our results so far. The eigenfunctions 
for the Alfven continuum are sine waves with an exponen- 
tial envelope. The wave length increases with 77 due to the 
increase in the Alfven speed. The eigenfunctions of the cusp 
continuum show a characteristic spike at the point where 
the field lines point radially outwards (and, less distinctly, 
at the equator for the closed field lines). This spike is a con- 
sequence of the fiattening out of the 77 profile of the partial 
hydrodynamic pressure. 



4.3 Eigenvalues 

We now evaluate the frequency spectrum for the magnetic 
mountain in Figs. [2] and |3] by computing lu'^ for polar (foot- 
point at ^ = 0.01 rad) and equatorial (footpoint at 6 — 1.6 
rad) field lines. Since B{R = i?,) decreases monotonically 
with ^, these two field lines bracket the continuous spectrum 
of the whole configuration. The results in section for the 
four field lines in Figs. [5] and [3] do not show any evidence for 
a spectrum that folds over onto itself, a feature i n some other 
MHD systems [cf. the discussion in ICocdbloed fc PoedtsI 
(|2004h ]. The spectrum folds onto itself when the charac- 
teristic speeds (Alfven and slow-magnetosonic in the plane 
parallel slab) have local extrema in the ^I" range considered. 

Fig. H displays lo'^ for Ma/M^ = 1.4, 1,0.6,0.1,0.01. We 
upscale all frequen cies to a realistic neutron s tar according 
to Lj'^ oc {ho/R,f (| Vigelius fc MelatosI [20083 ) . The Alfven 
(cusp) continuum is drawn in red (blue) , as a shaded interval 
with double-headed arrows, for node numbers n = 0, 10, 100. 
We stress that the gaps appearing in Fig. |3] stem from our 
choices of n; in fact, the whole spectrum lj^ ^ '^min is cov- 
ered without gaps by nodes n ^ 0. Furthermore, both the 
Alfven and the cusp continua for different n overlap. Overlap 
is common in some MHD systems, where multiple degenera- 
cies of the eigenfunctions occur. 

Let us compare the results in Fig. U] to the frequency 
spectrum of a plane parallel slab. We remind the reader that, 
in a one- dimensional system like the plane-parallel slab, the 
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71=0 n=1 n=30 



Figure 3. Field-aligned velocity eigonfunctions Jv^ of the cusp continuum for a magnetic mountain with accreted mass Ma = Mc = 
1.2 X IQ-^'Mq. The foot points of the field lines (D, ®, CD, and ® are at 6* = 0.1, 0.2, 0.6, 1 rad (left panel, tracing the field lines of the 
mountain in the meridional plane). The right panels graph Jv^ against arc length, rj, along the four field lines (top to bottom) for node 
numbers n = (left), 1 (middle), 30 (right), ho = 53.82 cm is the characteristic scale height (units of axes in left panel) and Cs = 10* 
cm s~^ is the sound speed, ri is normalized to the length of the field line, such that ^ ^ 1. 
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Figure 4. Alfven (red bands) and cusp (blue bands) continuous 
frequency spectrum (in Hz) for magnetic mountains with (from 
top to bottom) Ma/Mc = 1.4,1,0.6,0.1,0.01, computed for the 
node numbers (left bands to right bands) n = 0, 10, 100. 

spectra for difFerent n are defined by the dispersion rela- 
tions given in paragraph two of section |4]2] since the equi- 
librium values and thus the coefficients of equation (|28p do 
not depend on rj. The lower bound for tlie Alfven contin- 
uum is measured from Fig.|4]to be cjA,iow/(27r Hz) = 15.6 for 
Ma/Mc = 1. At tlie same time, the local Alfven frequency 
for the polar field line, lies in the range 3.39 uiA/i'iTT Hz) ^ 
55.38. Our aim is to pick the n = mode, so we assume 
ko ~ l/ho- In our two-dimensional system, and uj% are 
both functions of 77. The first grid cell is centered a.t = 0.01 
rad and hence determines the lowest polar field line that we 



can choose. This field line runs almost parallel to the z axis, 
i.e. r ~ const. Hence, the bottom row of (|28|) reduces to the 
dispersion relation for Alfven waves. Consequently, o^a.Iow 
lies in the range of the local Alfven frequencies. On the other 
hand, we can write down the dispersion relation for the cusp 
continuum assuming that A'' and B do not depend on 77, viz. 
uj^ = ujg — A/(NB^). The singular frequencies are offset by 
a constant due to the gravitational force. We compute uj{r]) 
for the polar field line to be 154 cj/(27rHz) 366. Again, 
we find that ci;c,iow/(27r Hz) = 189 (cusp continuum) lies in 
this range. 

For completeness, we display to for the four field lines 
in section 14.21 as a function of the node number n in Fig. [S] 
The Alfven (cusp) eigenfrequencies are plotted as plus (star) 
symbols. For field lines ® and ®, we find ujc > wa, consistent 
with the discussion in the previous paragraph. However, field 
lines (D and ® obey uja > i^c- The reason is that the latter 
fieldlines stay close to the surface. As a consequence, the 
constant offset, which is set by the gravitational force, is on 
average ~ 2.5 times lower for (D and ®, compared to ® and 
@. We find uja <^c for n — > for (D and ® , because A is 
neglible in this case and the eigenfrequencies depend on the 
node number as lo'^ oc with difFerent coefficients for uja 
and ujc- 

We are now in a position to explore the dependence of 
ojiow on the accreted mass Ma (Fig- Sjl and on the surface 
magnetic field B, (Fig. |6]). We do this by computing the 
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Figure 5. Alfvcn (star) and cusp (plus) cigcnfrcqucncics for a 
magnetic mountain with Ma = Mc and B* = 10^^ G for tlie four 
field lines ®— ® in section [4. 21 (four panels, labelled at top left) as 
a function of node number n < 10. 



Alfven and cusp frequency ranges employing the dispersion 
relations for constant field values established in the previous 
paragraphs and examine how they change with Ma and B, , 
respectively. The results are tabulated in table [T] 

The lower bounds for the Alfven (cusp) contin- 
uum are found to be [.tJA,iow/(27r Hz) — 191,15.6,13.7 
[t<JC,iow/(27rHz) = 73.7,189,183] for Ma/Mc = 0.01,1,1.4, 
respectively. It is interesting to compare Ma = O.OlAfc (early 
stage of magnetic burial) with Ma = Mc (middle stage). For 
Ma ~ O.OlAfc, the mass density is « 0.03 times the value 
at Ma = Mc and we find B{Ma = O.OIM^) ^ 2B{Ma = 
Mc), where B is averaged over field line ®. The gravita- 
tional term A is only moderately affected, with A{Ma = 

0. 01Mc) « 0.5A(Ma = Mc). On the other hand, we find 
B{Ma = lAMc) « 0.95B(Ma = Mc), N{Ma = lAMc) ~ 
1.24N{Ma = Mc), and A{Ma = lAMc) « lAA{Ma = Mc). 
This is the reason why ljiow peaks at M = Mc. For M = 

1. AMc, the magnetic field strength remains unchanged but 
the additional accreted matter increases A*'. 

We redo the analysis of the previous paragraph, 
this time varying the polar magnetic field _B* (Fig. 
[B]). The lower bounds for the Alfven (cusp) continuum 
are cjA,iow/(27r Hz) = 11.7,13.8,13.8 [a;c,iow/(27r Hz) = 
164,177,177] for B,/G = 10",10^^10" and Ma = Mc 
respectively. The continuous spectrum does not shift much 
as -B* varies over two decades. This is not surprising: Mc 
scales as Mc oc so that two mountains with a different 
-B* but the same Ma/Mc effectively share the same steady- 
state hydromagnetic structure and hence the same MHD 
spectrum. 



lO" 

//Hz 



Figure 6. Alfven (red bands) and cusp (blue bands) continu- 
ous frequency spectrum (in Hz) for magnetic mountains with 
(from top to bottom) polar magnetic field strength -B*/G = 
10^^, lO'^'^, 10^^, computed for the node numbers (left bands to 
right bands) n = 0, 10, 100. 



4.4 Rotational splitting 

One expects to find magnetically confined mountains on 
rapidly rotating neutron stars with 0.1 kHz < Q./2ti < 
0.7 kHz, e.g. in low-mass X-ray binaries, and accreting mil- 
lisecond pulsars. In such objects, the Coriolis force is an 
important factor in determining stability. We do not treat 
the Coriolis force in this article, because our main aim is to 
compare the analytically derived MHD spectrum with the 
output from zeus-mp simulations in the literature, which se t 
n = (|Pavne fc Melatosll2007l : I Vigelius fc Melatosll2008d). 
Howev er, the framework in section [3] and Hellsten fc Spiej 
l| 19791 ) can easily accomodate f2 7^ 0, as foreshadowed at the 
end of section 13.11 

Rotation splits the Alfven and cusp continua, possibly 
destabilizing them. We perform an order-of-magnitude cal- 
culation to estimate the influence of rapid rotation on sta- 
bility. Note first that the constant term A due to gravitation 
exceeds its centrifugal counterpart (yli,grav ~ GM^/R, 3> 
^i.rot ~ Jl^-R*), justifiying the ne glect of the centrifugal 
force in previous numerica l papers (jPavne fc Melatog|2007l : 
I Vigelius fc Melatos|[2"008cr ). For Q f 0, the two equations of 
couple via the Coriolis force (|Hellsten fc Spie^l 19791 ): 







and 



o-'^NB'^Jv'' 



i2anNr{B ■ Vr)v^ + 



B ■ V 



£2 



- jp 



-B ■ VJv 



e 



= cr^iVr^u"^ + \2aQ,Nr[B ■ Vr)Jv'^ 



(32) 



Vv^. (33) 



Clearly, the Coriolis force, which produces the second term 
in each equation, dominates the other terms. For (|32|l . we 
find \[2anNr{B ■ S/r)]/{a'^NB)\ - Q/a = 10^ > 1, and for 
the third term \ [2aQNr{B ■ Vr)]/[B ■ V{NB'^)(B^ + ■yp)B ■ 
VJ]| ~ (crn77^)/c? = 10*^ » 1. Equivalents, for dSHl), we 



n/a = 10^ > 1, and 



find \[2aD.Nr{B ■ VrJ)]/{a^Nr^) 

\[2aD.Nr{B ■ VrJ)]/{B ■ Vr^B ■ V)| ~ aQN/B^ = 10^ » 
1. A more general analysis including rotation is therefore 
needed in the future. 



5 DISCUSSION 

In this article, we compute semi-analytically the continu- 
ous part of the ideal-MHD frequency spectrum of axisym- 
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Table 1. Comparison of the lower bounds of the Alfven and cusp frequency spectra with the range computed under the assumption 
of constant field values (see text in section I4.3|l . The top and bottom panels record the eigenfrequcncies versus accreted mass Ma and 
surface magnetic field respectively. 
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Figure 7. Fourier transform t(f) of the mass cUipticity taken 
from ZEUS-MP simulations of magnetic mountains with Ma = Mc 
(solid), Ma = lAMc (dotted), and Ma = 0.6Mc (dashed). The 
continuous Alfven spectrum is clearly visible for Ma = Mc as a 
band starting at fa 11 Hz. We mark the theoretical analytic values 
of /a.Iow = 11-1, 16, 14 Hz for Ma/Ma = 0.6, 1, 1.4, respectively, 
as vertical lines. Note that the cusp continuum, with loq = 
177.07 Hz, lies outside the range of this plot, since it exceeds the 
Nyquist limiting frequency o^n = 44.5 Hz. 



metric, magnetically confined mountains. We find that the 
continuous spectrum covers all frequencies above a mini- 
mum cJniin- Furthermore, we find a;^;^ > in all the con- 
figurations (0.01 s; Ma/Ma 1.4) we study. This fur- 
ther substantiates two important properties deduced pre- 
viously from numerical simulations: (i) magnetic mountains 
are marginally stable, i.e., uj has zero imaginary part; and 
(ii) mountains relax hydromagnetically through the undu- 
lating submode of the thre e-dimensional Parker instability 
l|Vigehus fc Melatod [20083 ). which possesses discrete eigen- 
values only. We also find that, for Ma > Mc, when the mag- 
netic structure of the mountain substantially deviates from 
a dipole, oJmin for the Alfven and cusp continua depends 
weakly on Ma and B, . 

How do our analytic results compare with nu- 
merical simulations of oscillating magnetic mountains 



jPavne fc Melatosll2007l : IVigelius fc Melatosll2008j )? To an- 
swer this question, we perform an axisymmetric ideal-MHD 
simulation of a magnetic mountain with Ma ~ Mc, which 
is perturbed slightly at t = 0. The simulation is performed 
using the parallel ideal-MHD solver zeus-mp (|Haves et al.l 
I2OO6I ). Fig. [7] shows the Fourier transform e(/) of the mass 
ellipticity e{t), which is proportional to t he mass quadrupole 
momen t of the mountain [see Eq. (2) in lVigelius fc MelatosI 
(|2008j )]. We choose e to compute the spectrum as (i) it 
is directly measurable from future gravitational wave data 
(|Vigehus fc Melatod[2008ah and (ii) it is an integrated quan- 
tity sampling A'' and B everywhere, so it is sensitive to all 
global oscillation modes. 

Besides the constant offset at / = 0, the simulation 
with Ma = Mc (solid curve in Fig. [7| exhibits a continuous 
spectrum with a lower boundary at / ~ 11 Hz. This is con- 
sistent with the analytic theory in section 2] which yields 
/a, low ~ 16 Hz, indicated as a vertical line in Fig. [T] The 
width of a frequency bin is 5f = 1.45 Hz, so we conclude 
that the lower boundary of the simulated spectrum almost 
coincides with /a, low It appears that the lowest modes of 
the continuous Alfven spectrum are indeed excited in ZEUS- 
MP simulations, substantiating the claim that global MHD 
mountain oscillations are Alfvenic (|Pavne fc Melatosll2007l ). 
For Ma = lAMc (dotted curve in Fig. [7|, we see a continu- 
ous spectrum at / > 7 Hz, albeit less distinctly than in the 
Ma = Mc case. The Ma = O.QMc simulation (dashed) shows 
a peak at / = 23 Hz and no obvious continuous spectrum. 

We do not attempt a unambiguous identification of the 
oscillation modes seen in the above simulation. As an in- 
tegrated quantity, the behaviour of e is determined by a 
global superposition of different eigenmodes, both, continu- 
ous and discrete. These eigenmodes ar e stochastically ex- 
cited through numerical inaccuracies (|Vigel ius fc MelatosI 
l2008d ) in the simulation displayed in Fig. [7] and identifying 
the frequency spectrum of e with the underlying eigenfunc- 
tions is hence difficult. The calculations undertaken in this 
article are but a first step towards a thorough investigation 
of the complete eigenvalue problem. 

In the context of future gravitational wave observa- 
tions of magnetic mou ntains in low-mass X-ray binaries 
(|Melatos fc PavndbOOSl ) . the MHD oscillation spectrum en- 
ters the gravitational wave signal through sidebands and 
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broadening near the Fourier peaks at / « and 2/«, where /<, 
is the s pin frequency of the neutron star. lVigelius fc Metatod 
(|2008al ) showed that the sidebands can be observed in princi- 
ple with next-generation interferometers. However, in order 
to exploit these features fully to probe the physics of surface 
magnetic fields on accreting neutron stars, we need to know 
more about how the oscillations are excited (e.g. by the vari- 
able accretion torque) and damped. The analytic technique 
in this article is a useful tool for such investigations. 

The oscillation modes may be perpetually re-exci ted, 
e.g. through starquakes, variable accretion torques (jLail 
I1999I). and possibly cycl o nic flows during type I X-ray bursts 
(jSpitkovskv et al.ll20o3 ). IVigelius fc MelatosI (|2008al l show 



that a perturbation of the fluid density causes a fractional 
change in the signal-to-noise ratio d of 5d/d ~ n/N. Sim- 
ilarly, in a (very crude) model where the pulse shape is 
determined by the positions of the footpoints of the mag- 
netic field lines at the polar cap boundary, the fractional 
change in pulse parameters (e.g. full-width half-maximum) is 
comparable to the fractional perturbation amplitude. Unfor- 
tunately, the excitation mechanisms are poorly understood 
and it is unclear if the amplitude of the perturbation is suf- 
ficient to cause observable features in the gravitational wave 
spectrum. Ultimately, gravitational wave observations will 
yield valuable information about the underlying excitation 
phy sics. 

iGlampedakis et all ()2007h investigated ideal-MHD 
modes in magnetars, taking into account the coupling be- 
tween the fluid core and the elastic crust. They found that 
global core-crust modes can explain quasi-periodic oscilla- 
tions (QPOs) observed during giant flares in the soft gamma- 
ray repeater s SGR 1806-20 and SGR 1900-1-14. In contrast, 
iLevinI (j2006h argued that continuous coupled modes in mag- 
netars decay too rapidly to account fo r the observed QPOs. 
This debate was reviewed recentlv bv I Watts &: Strohmaveil 
l|2007l ). Using a series expansion, iLed (|2007l . l2008h calcu- 
lated the discrete eigenmodes of magnetars. While low fre- 
quency QPOs can be identified with fundamental toroidal 
torsional modes, higher frequency (100 < / < 1000 Hz) can 
be attributed to a variety of modes, such as spheroidal shear 
modes or core/crust interfacial modes. 

The aim of this paper is to interprete analytically 
and physically the magnetic moun tain oscillations seen in 
nonrotating nu merical simulations (jPavne fc Melato3l2007l : 
IVigelius fc Melatosi.2008ai 'l . For rapidly rotating objects with 
f2/27r > 0.1 kHz, like accreting millisecond pulsars, the con- 
tinuous spectrum is strongly modified by rotation, as shown 
in section We will calculate the rotational splitting in a 
forthcoming paper. 
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